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Abstract 

In this note we consider several kind of partition functions of one- 
dimensional models with nearest - neighbor interactions /„, n e Z and 
spin values ±1. We derive systems of recursive equations for each kind 
of such functions. These systems depend on parameters /„,n G Z. 
Under some conditions on the parameters we describe solutions of the 
systems of recursive equations. 
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1 Introduction 

The properties of physical, biological and many other systems can be de- 
scribed by differential and recursive equations; the latter are also called dis- 
crete dynamical system (see e.g. |2], [3], [IQ]). Also systems of non-linear, 
higher dimensional recursive equations arise in solving many different prob- 
lems (see e.g. [T], [S], [B], [2], [II], [I2])- But theory of the systems of recursive 



equations is not developed enough. So for each such a system one has to 
use a specific argument which is suitable for soling the system. 
In the paper we consider the Hamiltonian (energy) 

-^(^) = X! Ix'i-aix-l)jtaix)-, (1) 

l={x-l,x):x€Z 

where Z = {...,—2,-1,0,1,2,...}, a is a function (configuration), a : Z ^ 
{ — 1,1}, (the set of all configurations a is denoted by $7 = {—1,1}^) and 
Ix ^ R for any x G Z. 

In statistical physics the Hamiltonian ([1]) is called an one-dimensional 
(ID) model. 

Let us consider a sequence = [— n, n], n = 0, 1, ... and denote AJ^ = 
Z\ A„. Consider a boundary condition o"!^^ = o"ac = {cj{x) = +1 : x G AJ^}. 
The energy H^{a) of the configuration a in the presence of the boundary 
condition cri^^ is expressed by the formula 

(2) 

l={x — l,x):x(I^A„ 

The Gibbs measure on 0^ = {—1, 1}^" with respect to the boundary 
condition cri^^ is defined by 

//+^(<t) = Z-\n,P,+)exp{-PH+ia)), (3) 

where P = T^^, T > 0— temperature and Z(n, /?,+) is the normalizing 
factor (statistical sum or partition function): 

Zin,P,+)= J2 expi-PH+iip)). (4) 

Note that the probability (with respect to measure /i^^) of a subset Q'^ of 
Cln is defined by 
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where Z'(n,P,+) is called a "crystal" partition function: 

Z'{n,f3,+)= J2 exp(-/?if+(V)). (6) 

So to define the Gibbs measure and probability of an event of the system 
one has to compute the partition functions. If /u^ = lim„^oo IJ-np exists then 
it is called a limit Gibbs measure. A limit Gibbs measure for a given type 
of interaction (energy) may fail to be unique this means that the physical 
system with this interaction can take several distinct equilibria i.e there is 
phase transition. 

Note that (see |4j, p. 95) for the model ([I]) on TV = {1,2,...} it was 
shown that there occurs a phase transition iff J2n>i < oo. In [8j using 

a contour argument it has been proven that for that model ([1]) the phase 
transition occurs if + In+k > ^ for any n G Z, k £ N. 

In this paper we consider some (crystal) partition functions of the model 
and give the system of recursive equations for the functions. Under some 
conditions on parameters of the model we describe their solutions. 

2 Partition function of "-|-" and "±" -boundary 
conditions 

Consider two type of partition functions: 

Z+= exp{-PH+{an)}, (7) 

Zt= exp{-/3i?±((7„)}, (8) 

cr Yi Gf^n 

where is defined by ([2]) and 

H^ian) = H+ian) + I-n^{-n). (9) 
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In this paper for the simplicity assume 

In = I-n+i, for any n € Z. (10) 



Proposition 2.1. // condition \10\) is satisfied then the partition functions 
^ and |3j have the form 



ill) 



Proof. Under the condition (jlOp we get = and Z^ = Z^ . Now from 
([7|) , ([8]) we obtain the fohowing system of recursive equations 



Z+ = (1 + e-2/^^"+0^+-i + 2e-'3^"+^Z±_i, 
Z± = (1 + e-2/3/.+i)Z±_, + 2e-/^^"+^Z+_, 

Putting Xn = Z^ — Z^ and = Z^ + Z^ from (fT2]) we get 



(12) 



(13) 



X„ = (l-e-'3^"+i)2x„_i, 

Yn = (l + e-^^^+O'^n-l. 

The equaUties Xq = Z+ - = (1 - e-^^i)2, = (1 + e-^^i)^ with 
(fT3D imply 

n n 

X„ = - e-^^»+0^ = 11(1 + e-^^'+i)2. 

4 = 1=0 

Hence we get pT]) . □ 

For example, in a case of the usual Ising model i.e. Im = I, Vm G Z 
from pT]) denoting r = exp(— /3I) we get 

Z+ = i ('(I + r)2{"+i) + (1 - r)2("+i) ) , 
Z± =!((! + r)2("+i)-(l-r)2(«+i) 



Using these equalities (for usual Ising model) it is easy to see that 

— >■ 1, if n ^ oo. 

This means that for the Ising model the partition functions and Z^ are 
asymptotically equal. This gives in fact uniqueness of limit Gibbs measure 
for the ID Ising model. Such an asymptotical equality is true if Im is a 
periodic function of TTi i.6 Irff^j^p — I^a for some p ^ 1 Siiici all tyi G A^. 

3 Crystal partition functions 

In this section we are going to describe the crystal partition functions. 

Denote ^rn,n = { — 1,1}'™'"^, where [m,n] = {m,m + l,...,n}, m,n € 
Z, n > m. Put 

Ne{a) = |{x G [m, n] : a{x) = e}\, e = ±1, 

where IS*! is the cardinal of the set S. For r = 0, 1, ...,n — m + 1 consider 
the following crystal partition functions: 

Z^m:n= E ^-'""^^"^^ ^ = -'+ (14) 

o-Gf2m,„:Af+((T)=r 



Note that Z^f^ = Z+f;^. Denoting = Z^f^ and y;;„ = Z^f^ 

from (jl4p and (jlSp one easily gets the following system of (multi-variable) 
recursive equations 

(16) 

yr _ yr i -/37„+i V^-l 

m,n m,n—l ' ^ m,n— 1' 

where r = 0, 1, n — m + 1, m, n G Z, n > m. 



The system (fT6|) can be reduced to a recursive equation with respect to 
X^ yj. Indeed, from the first equation of (fTUj) we get 

^ m,n—l ^ \^m,n m,n—lj ' \^') 

Now from the second equation of p6]) using pT|) we get 

^^m,n ^^m,n—l ' ^ \^^m,n—l ^^m,n—2i ' ^ m,n— 2' 

(18) 

where r = 1,2, n — m + 1, n > m and 

vO _ 1 _ p-/3(/m+-fm+i) ;v-o _ 1 

_ p-Z^C/m+^fm + l) _l_p-/3(-/'m + l+-fm + 2) Tt^2 _ „-/3(/m+-fm + 2) rnP7, 



Iterating (jlSl) one can obtain an expression for Then using (11^ 

one can find Y^^. But these expressions would be in a very bulky form. 

Now we shah illustrate such an expression for the Ising model i.e. = 
I, m G Z. In this case the recurrence equation (|18p becomes more simple 



x: = x:z[ + x:_, + u- (19) 

where X^ = X^ „^ with ni — m = n, x = e~^^^. 

It is easy to see that for r = 0, 1, 2, 3, 4 the solutions are 



X0 = 1, Xi, = nx, Xl = {n-l)x + 



(n - 2)(n - 1) 
2P 
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Xf, = (n - 2),, + („ - 2)(„ - 3),= + (~-2)(»-3)(„-4) ^3_ 
A- = (,. - 3), + ^"' -^f-^' x' + (n-3)(„-4)(,.-5) ^3^ 

(n-3)(ra-4)(n-5)(n-6) 4 
4! ^ ' 
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here we used the foUowing formulas 

n 2^ 1 

J2f = ^(2n=^ + 3n2 + n), ^ j^^ = -(n^ + + n^). 
i=i i=i 

Note that has a form 

K = alnX + a2,nX^ + ... + <,nX''. 

For the coefficients ^, < k < r < n, using (jl9p we obtain the following 
system of recursive equations 

f 7^ I T — 1 T — 1 

"l,n — ^l,ra-l + '^l,n-l ~ ^l,n-2' 

afe,n = «fc,n-l + + " arn-2' ^ = 2, 3, r - 1, (20) 

r r I 1 

The following lemma gives solution to (|20p 

Lemma 3.1. Solution of the system of recursive equations (2D() is 

< k < r < n. (21) 



''k.n 



n — r + 1\ Ir — 1 
k U-1, 



Proof. We shall use mathematical induction (cf. with [7J pages 148-150). 
Let Am denote all cases of ()2ip with n + k + r = m. Formulas given above 
for X^, r = 0, 1, 2, 3, 4 show that the formula (|2ip is true for small values of 
m. Assuming that Am holds, we are to prove Am+i that is, equation ([20]) . 
for any integers n, r and k whose sum is m + 1. Since RHS of equation (j20p 
contains terms with n + r + k < m using the assumption of the induction 
for each term of RHS of ([20]) we get dH]). □ 



Thus the solution of ([T9|) is given by 



K=tr:^%-_%^ (22) 



k=l 
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Remark 1. Note that for X = 1 ('i-^- there is no interaction) the solution 
of < f7P|) is = (") . Using for x = ^ we obtain the following property 
of binomial coefficients 



This identity is known as the convolution identity of Vandermonde. 

Since interaction (parameter /) of the ID Ising model is translation- 
invariant (does not depend on the points of Z), the unknown functions 
-^mm ^mn t^c systcm (fT6]) depend on n — m and r only (see (fT9|) ). 
consequently, instate of n — m we can write n. Summarizing the results for 
the Ising model we have 

Theorem 3.1. For the Ising model the solution of the system of recursive 
equations U6\) is 



where n stands for n — m. 
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